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We present results of theoretical description and numerical calculation of the dynamics of molec- 
ular liquids based on the Reference Interaction Site Model / Mode-Coupling Theory. They 
include the temperature-pressure(density) dependence of the translational diffusion coefficients 
and orientational relaxation times for acetonitrile and methanol in water at infinite dilution. 
Anomalous behavior, i.e. the increase in mobility with density, is observed for the orientational 
relaxation time of methanol, while acetonitrile does not show any deviations from the usual. 
This effect is in qualitative agreement with the recent data of MD simulation and with experi- 
mental measurements, which tells us that presented theory is a good candidate to explain such 
kind of anomalies from the microscopical point of view and with the connection to the structure 
of the molecules. 

Introduction 

Dynamics in solutions has been a central issue in the physical chemistry. It has been 
so because liquid dynamics is an important factor which determines the rate of chemical 
processes in solution including chemical reactions and protein folding. A variety of ex- 
perimental methods have been devised to observe liquid dynamics, which now cover the 
dynamic range from femtoseconds to days. 

Theoretical studies of liquid dynamics and relaxation have been dominated during the 
past century by the two continuum models which had been established essentially in the 
18th century: the electrodynamics and the hydrodynamics [1]. The monumental achieve- 
ments made in the field of liquid dynamics by the great scientists such as A. Einstein [2], 
P. Debye [3], L. Onsager [4,5] and M. Born [6] are essentially based on the continuum 
theories: the dielectric relaxation, the ion mobility and its concentration dependence, 
and so on. Such theories represented by the Stokes-Einstein-Debye (SED) law rely on the 
boundary conditions as well as phenomenological use of macroscopic constants such as the 
viscosity and dielectric constants in order to realize "chemistry" specific of the problem 
in concern. Typical examples of the boundary conditions are the "slip" and "stick" ones 
often employed in equations based on the hydrodynamic theory. 
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The continuum models have played crucial roles in describing physics taking place in 
solution with the spatial as well as temporal resolutions which are low enough so that 
molecular details of the process can be neglected. There have been gaps between the 
phenomenological model and molecular processes. The gaps have been filled by using 
such adjectives as "effective". For example, the conductivity of a simple ion in aqueous 
solutions at infinite dilution depends on their size in entirely opposite manner to what is 
predicted from the Stokes law, if one employs the crystal radii as the ion size. In order to 
fill the gap between the phenomenological model and the molecular process, people have 
invented an "effective" radius or the Stokes radius, and interpreted the phenomenon in 
terms of the enhanced Stokes radii due to "hydration" . So called "fractional exponent" , 
e.g. [7], often used in the interpretation of the viscosity dependence of a relaxation rate 
in liquids is another example which tries to compromise the contradiction between phe- 
nomenological models and experiments. If an experimentalist tries to plot a relaxation 
rate against viscosity, he will often face serious violation of a law predicted by a phe- 
nomenological theory. It has been a common maneuver in such cases to use "fractional" 
dependence on viscosity. However, the use of "fractional" dependence does not add any 
information to describe molecular process actually taking place in solutions, but just hides 
the breakdown of the phenomenological theory. Recent experimental techniques, which 
are dramatically improved both in resolution and sensitivity, have demonstrated unam- 
biguously breakdown of the phenomenological model. Let us refer to two contributions 
which exhibit such breakdown. Nakahara [8] and co-workers have shown that the rota- 
tional dynamics of a solute in a variety of solvents does not correlate with the viscosity 
of the solvent but rather with specific interactions between solute and solvent. Terazima 
and his coworkers [9] have found that the diffusion constant of radical species is roughly a 
half of that of their parent molecules, which are about the same size. We suppose no more 
words are necessary to convince people what may happen if one tries to describe more 
complicated dynamics in solution, e.g., conformational change of protein and chemical 
reactions, in terms of the phenomenological models. 

Are there, then, any hope for theories of solution dynamics to break through the old 
regime established by the great scientists? We say "yes" , if one relies on the two theories in 
the statistical mechanics developed in the last century: the generalized Langevin equation 
(GLE) [10] and the RISM theory of molecular liquids [11-13]. The generalized Langevin 
theory describes the time evolution of few dynamic variables in concern as a function 
of the representative point in the phase space. All other variables which are not under 
explicit concern are projected onto the dynamic variables of interest with the help of 
a projection operator. The projection leads to an equation which looks similar to the 
phenomenological Langevin equation containing the frictional force proportional to the 
rate of change of the dynamic variables, and the random force, which are related to each 
other by a fluctuation-dissipation theorem. If one chooses as the dynamic variables the 
density and the conjugated flux of sites or atoms of molecules in liquids, the theory gives 
an equation for the dynamic structure factor of atoms, which describes the time evolution 
of the site-site density pair correlation functions [14]. Results from the RISM theory, the 
site-site static structure factor and the direct correlation functions, are naturally plugged 
into GLE in order to evaluate not only the initial condition of the dynamics or the static 
structure factor, but also the frictional force as well as the collective frequency which 
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concerns the frequency of the site-site density fluctuation. A crucial development of the 
theory is rather conceptual, not mathematical, in the sense that it has provided a new 
concept to view dynamics of a molecule in solution, which is quite different from the model 
traditionally exploited in the field. The new model sees the molecular motion in liquid as a 
correlated translational motion of atoms: if two atoms in a diatomic molecule are moving 
in the same direction, then the molecule as a whole is making a translational motion, 
while the molecule should be rotating if its atoms are moving in opposite directions. The 
view is different from the rotational-translational model traditionally developed based on 
the angular coordinates [15-17]. 

The new theory of liquid dynamics has been successfully applied to a variety of relax- 
ation process in solution including the collective excitations in water [18], ion dynamics in 
dipolar liquids [19], dynamical Stokes shift [20,21], pressure dependence of the transport 
coefficients [22] as well as dielectric relaxation spectrum in water [23], and so forth. In 
the present proceeding, we report our latest study on dynamics of a molecule in solution 
at infinite dilution and its temperature and density dependence. 



Theory 

Equilibrium structure 

In this work, we use the dielectrically-consistent reference interaction-site model (or 
DRISM for brevity) formalism [24] for the system. The main equation here is the site-site 
Ornstein-Zernike equation (SSOZ) written as 

ph(k)p-x(k) = [&(k) + x{k)]c{k)[u(Q (1) 

where p is the diagonal matrix of number density of molecular species, and h(k), c(k) 
and u)(k) are the total, direct and intramolecular correlation matrices, respectively, in 
the reciprocal space. The matrix is determined by the dielectric properties of the 
system. 

In the case of infinite dilution limit, i.e. when concentration of one species (called a 
solute) tends to zero and concentration of other (called a solvent) essentially determines 
the total density, the DRISM equation can be decomposed as follows: 

h vv (A;) = [w v (A;) + D v (A;)p v ]c vv (A:)[w v (A;) + p v h w (Jfe)] + D V (A;), (2a) 
h uv {k)=w u {k)c uv {k)[w v {k) +p v h vv (fc)], (2b) 
h uu (A;)=w u (A;)c uu (£;)w u (£;) + w u (A;)c uv (A;)p v h vu (A;), (2c) 

where w(k) = u>(fe)p -1 , £)(&) = p _1 %(A;)p~ 1 , and superscripts "u" and "v" refer to 
"solute" and "solvent" , respectively. Equations (J2J) are solved with the hypernetted chain 
(HNC) type of closure specified as 



g ai [r) = exp 



retry \y ) 7 / \ / \ 
+ h a Jr) - c a Jr) 



(3) 

where functions without tilde are those in the real space, a and 7 are site labels, and 
g aj (r) = [hary{r) + 1] is the site-site radial distribution function. Interaction potential 
between sites a and 7 is denoted as Q7 (r), (3 = 1/k-eT with &b and T being the Boltzmann 
constant and the absolute temperature, respectively. 
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The site-site mode-coupling theory 

The site-site intermediate scattering function of neat solvent and its self-parts are obtained 
by the generalized Langevin equation / modified mode-coupling theory for molecular liq- 
uids as is described in the previous study by Yamaguchi et al. [22] . Generalized Langevin 
equations for the neat solvent are given by 

F(k,t) + k 2 J(k)-S- 1 (k)-F(k,t) + f dr K(k, t - r) ■ F(k, r)=0, (4a) 

Jo 

F s (k,t) + k 2 J(k)-w-\k)-F s (k,t)+ [ dTK s (k,t-T)-F s (k,r)=0. (4b) 

Jo 

Here, the site-site intermediate scattering function and its self-part in the time domain, 
denoted as F(k,t) and F s (/c,t), respectively, are given by 

F ai (k, t)=^(p* a (k, t = 0) Pl (k, t)), (5a) 

K r (k,t)=±(pZ(k,t = 0)f h (k,t))', (5b) 

where N is the number of solvent molecules, p a (k) is the density field of a-site in the 
reciprocal space, and (• • ■) means the statistical average, the suffix "s" means that cor- 
relations between the quantities of different molecules are neglected. Connection with 
equilibrium properties is expressed in the relation between the site-site static structure 
factor S(k) = F(k,t = 0) and direct and intramolecular correlation functions, which is 

S ai (k)=±(p* a (k,t = 0K(M = 0)) = [w v (A0 +ph vv (k)} . (6) 



J a-f 



The site-current correlation matrix J(k) is defined in a similar way as 

J aj (k) = ±(j* !Z (k,t = 0)j„(k,t = 0)), (7) 

where z-axis is taken to be parallel to the k vector. The mode-coupling expressions of 
the memory function matrices, denoted as K.(k,t) and K s (/c, t), are given by 

[j- 1 (A;)-K MC T(A;,t)] a7 =^p j dq {g, 2 L6(g)-F(g,t)-c(g)J a7 F a7 (|k-q|,t) 

- Qz(k - q z ) [c{q) ■ F(q, t)\ ai [F(|k - q|, t) ■ c(|k - q|)J a7 } , (8a) 
[J- 1 ^) • K S MCT (M)W=7^ j d q g 2 L6(g)-F(g,t).c(g)j a7 F^(|k- q |,t). (8b) 

According to the recipe by Yamaguchi and Hirata [25] , memory functions for the self-part 
are given by the linear combination of the corresponding mode-coupling memory functions 

as 



mk, t) . j(k)] ai = J2 (^L z z m2 zz m A^ ? - r '- r1 ^ 

mi t 2,3&{x,y,z} 

x [J-^^.K^^O]^, (9) 
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where Z^ lTn2 and u Z m stand for the orientation-dependent site-site velocity correlation 
matrix and the unitary matrix that describes the rotation between the molecular- and 
laboratory-fixed coordinates of molecule i, respectively [25]. The collective-part of the 
memory function (neglecting the orientational correlation between different molecules) is 
given by 

K(k, t) = K MCT (k, t) + K S (A;, t) - K S MCT (A;, t). (10) 

The time-evolution of the self-part of the solute-solute site-site intermediate scattering 
function is governed by the equations similar to Eqs. (|4b|) . (j8bj) and (jHJ) as 

F uu ' s (k,t) + k 2 3 u (k) ■w u - 1 (fc) -F uu ' s (k,t) + [ drK uu ' s (fc,t-r)T uu ' s (fc,r) = 0, (11) 

Jo 

[J™>-\k) ■ KZ&(k,t)U = j£y 3 j dqq 2 z [c uv ( q ) ■ F vv (q,t) ■ £ vu (g)j a7 F-/(|k - q|, t), 

(12) 

[K^(k, t) • J uu (fc)] Q7 = ^ZZ^ZruA^'^-^) 

^1,2,3 &{x,y,z} 

x [j^-\k)-K£&(k,t)]. (13) 



Diffusion coefficient and reorientational relaxation time 

Based on the Green-Kubo formula, the translational diffusion coefficient D is obtained 
as [26, 27] 

1 f°° 

D = -J q dt Z aj (t), (14) 

where Z(t) is the site-site velocity autocorrelation function, which is described in terms 
of the self-part of the intermediate scattering function F s (k,t) as 

z «r(*) = ^Z)<Va(0)-v 7 (t)>' = -jim ^F^(k,t). (15) 

i 

In Eq. ()14j) . a and 7 refer to any two sites in a molecule which are bound by chemical 
bonds. The expression for the single-particle reorientational time is also described in terms 
of the site-site velocity autocorrelation function [22,27,28]. In this work, we restrict our 
discussion to the rank-1 reorientational relaxation of the dipole moment /j, given by 

Mi = J> Q if, (16) 

a 

where z a is the charge of the site a. The first-rank reorientational correlation function 
C M (i) is defined as 

c * {t) = mm • (17) 
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Substituting Eq. (fTBj) into Eq. ([T7]). C M (t) is related to the site-site velocity autocorrela- 
tion function Z(t) as 

The reorient at ional correlation time of the 1st rank, r, is then defined as the time- 
integration of Cpit). 



Description of the models 

We performed explicit calculations for two popular systems, i.e. acetonitrile (CH3CN) in 
water and methanol (CH3OH) in water, both in the case of infinite dilution. As for the 
structure and the intermolecular potential of water we employed a model of the extended 
simple point charge (SPC/E) [29]. We also put the Lennard- Jones (LJ) core on the 
hydrogen atoms in order to avoid the undesired divergence of the solution of the RISM 
integral equation. The LJ parameters of the hydrogen atom, the depth of the well and 
the diameter, are chosen to be 0.046 kcal/mol and 0.7 A, respectively. 

In acetonitrile and methanol the methyl group is considered to be a single interaction 
site (Me) located on the methyl carbon atom. So that both chemical compounds consist 
of three sites interacting through the pair potential [30, 31] 

<j>{Ti,rj) = <f>(rij) = \ 4e «/3 

a,/3 { 

i.e., LJ plus Coulomb. Here a and (3 label sites on molecules % and j; r ia ^p = \r ia — r^l; 
parameters e a p and a a/ 3 are LJ well-depths and LJ diameters defined as e a p = ^Jt a e$ and 
&af3 = (cq + cr / a)/2, respectively, with a a being the LJ diameter of a single site. Point 
charges for acetonitrile were chosen to reproduce electrostatic potential obtained in an 
ab initio calculations [30]. Numerical values of parameters of the site-site interaction 
potential (flTj|) and masses of sites are specified in Table [T] 



cr. 
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o-, 



a/3 



+ 



Tia,jf3 



(19) 



Table 1: Parameters of the site-site interaction potential (|19|) : mass and charge are in atomic units, ctlj 
in A, and £lj in 10~ 14 erg/molec. 
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16.0 


1.008 


1.008 


15.024 


12.0 


14.0 


15.024 


16.0 


1.008 


z 


-0.8476 


0.4238 


0.4238 


0.269 


0.129 


-0.398 


0.265 


-0.7 


0.435 


CLJ 


3.16 


0.7 


0.7 


3.6 


3.4 


3.3 


3.74 


3.03 


1.0 


£LJ 


1.084 


0.3196 


0.3196 


2.64 


0.6878 


0.6878 


1.4525 


1.1943 


0.3196 






a Ref. [29] 






b Ref. [30] 






c Ref. [31] 





Information about bond length can be deduced from Cartesian coordinates (x, y, z) of 
sites indicated in Table |2] (principal axes and the origin can be taken arbitrarily). In 
calculations, the temperature of the system is varied from 258.15 to 373.15 K, and the 
number-density from 0.03008 to 0.04011 molecules/A 3 , where the number-density of water 
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Table 2: Cartesian coordinates of sites. 

water acetonitrile methanol 

site (x,y,z), A site (x,y,z), A site (x,y,z), A 

O (0, ,-0.0646) Me (0, 0, 1.46) Me (-1.4246, , 0) 

Hi (0, 0.8165, 0.5127) C (0, 0, ) O(0 ,0 ,0) 

H 2 (0,-0.8165, 0.5127) N (0, 0,-1.17) H ( 0.3004, 0.8961, 0) 



at the ambient conditions is 0.03334 molecules/ A 3 . Connection of these parameters with 
thermodynamic pressure is shown in Table El (except for the metastable regions where 
we do not have reliable data). Temperature/density dependent dielectric permittivity e 

Table 3: Density-pressure correspondence for water according to [32]. Pressure is given in MPa. 

p, g/cm 3 n, A" 3 273.15 K 298.15 K 373.15 K 

(19 0.03008 = = = 

1.0 0.03334 0.4085 6.6914 100.64 

1.1 0.03676 257.20 296.20 451.51 

1.2 0.04011 689.03 760.76 993.38 



for water used in numerical calculations has been evaluated as a physical solution of an 
empirical nonlinear equation presented in [33]: 

If 1\ if 9.32 ■ 10 4 (1 + -MA \ , s 

e-- l + -)=-17+ v 9 , 20 

where v is a molar volume in units of cm 3 /mol, and T is thermodynamic temperature in K. 
This equation was also used in such temperature/density points where no experimental 
data exist. 



Numerical methods 

From the generalized Langevin equation / modified mode-coupling theory and the DRISM 
/HNC integral equation theory, the diffusion coefficients and the reorient at ional relaxation 
times of solute molecules in solution can be obtained based solely on the information about 
molecular shapes, inertia parameters, intermolecular interaction potentials, temperature 
and density. First, we calculated the site-site static structure factor by the DRISM/HNC 
equation using the intermolecular interaction, molecular shape, temperature and density. 
In order to improve the convergence of the RISM calculation, we used the method of the 
modified direct inversion in an iterative space proposed by Kovalenko et al. [34]. From the 
static site-site structure factor, we calculated the site-site intermediate scattering function 
using the site-site generalized Langevin equation / modified mode-coupling theory. The 
generalized Langevin equation is time-integrated numerically. The time-development of 
the correlation functions in the k — > limit is separately treated by the analytical limiting 
procedure of the theoretical expressions. In the numerical procedure, the reciprocal space 
is linearly discretized as k = (n+ |)AA;, where n is the integer from to iV fe — 1. Values of 
Ak and Nk are 0.061 A -1 and 512, respectively. The diffusion coefficient D was calculated 
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from the asymptotic slope of the time dependence of the mean square displacement and 
the orientational relaxation time r was analyzed using the rotational autocorrelation 
functions. 



Results and discussion 

Figure ^ shows the density dependence of normalized diffusion coefficients D/D Q , with D 
being the diffusion coefficient at the lowest density, for both acetonitrile (left part) and 
methanol (right part) at four temperatures. Corresponding results for the orientational 
relaxation time r are shown in figure El It has been found that diffusion coefficients 
of both acetonitrile and methanol at these temperatures decrease monotonically with 
increase of density except for methanol at the lowest temperature, where one observes a 
typical anomalous behavior: first diffusion coefficient increases, then decreases again. At 
the same time, behavior of orientational relaxation time is different for acetonitrile and 
methanol at all studied temperatures, except at the highest. One can easily observe that 
r for acetonitrile exhibits the monotonic increase with density and monotonic decrease 
of its absolute value with temperature [figure 12(a)], while r for methanol at the same 
temperatures first decreases with density and then increases [figure 12(c)] . Its absolute 
value is also getting smaller with temperature as in the case of acetonitrile. In such a 
way, the density dependent orientational relaxation time for dissolved in water methanol 
has an absolute minimum somewhere in between 0.9 g/cm 3 and 1.1 g/cm 3 for water and 
the depth of this minimum is seen more clearly at lower temperatures. These tendencies 
are observed sharply in figures 12(b) and [2(d) for acetonitrile and methanol, respectively, 
where we plot normalized values of r as t/tq with To being the orientational relaxation 
time at the lowest density. Obtained results are in qualitative agreement with the outco- 





p, g/cm 

Figure 1: Density dependence of diffusion coefficients for acetonitrile (■) and methanol (A) in water at 
different temperatures. 

me of the MD simulation by Chowdhuri and Chandra [35] performed for similar systems 
at 258 K and 298 K, and with experimental measurements by Wakai and Nakahara [8] 
who observed likewise behavior for several different systems, including acetonitrile, at 
room temperatures. 

The difference in behavior of acetonitrile and methanol can be made clear from the 
point of view of their molecular structure. To begin with, let us remind that the anomalous 
density dependence of the molecular mobility in water has been explained by Yamaguchi 
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Figure 2: Density dependence of the oricntational relaxation times for acetonitrile (■) and methanol 
(A) in water at different temperatures. 



et.al. [22] based on the facts of an almost spherical repulsive core for water and the strong 
short-range intermolecular Coulomb interaction called "hydrogen bonding". Following 
the idea and owing to the models used in present manuscript, one can conclude that since 
acetonitrile cannot make the hydrogen bonding and methanol does, this fact may be a 
good candidate to serve as a reason or cause of justification of the anomalous density 
behavior of methanol (or possibly any protic solution) and usual one of acetonitrile (or 
possibly any aprotic solution). 

Summary 

In present work we have calculated the temperature/density dependence of the transla- 
tional diffusion coefficients and rank-1 orientational relaxation times for acetonitrile and 
methanol in water at infinite dilution using the site-site generalized Langevin equation / 
modified mode-coupling theory and the DRISM/HNC theory. Calculations show anoma- 
lous density dependence of the orientational relaxation time for methanol in water which 
is consistent with the results of experimental observation and MD simulation. On the 
other hand, similar computation for acetonitrile does not exhibit deviation from the usual 
behavior which is also in agreement with experiment and MD simulation. Such a differ- 
ence is explained based on the molecular structure of solvents, in particular their ability 
(for methanol) and inability (for acetonitrile) to create hydrogen bonds. More studies will 
be required for better understanding of the molecular dynamics in such systems. 
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